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I discuss some non-perturbative aspects of hot gauge theories as related to the 
unscreened static magnetic interactions. I first review some of the infrared diver- 
gences which cause the breakdown of the perturbation theory. Then I show that 
kinetic theory, as derived from quantum field theory, is a powerful tool to construct 
effective theories for the soft modes, which then can be treated non-perturbatively. 
The effective theory at the scale gT follows from a collisionless kinetic equation, of 
the Vlasov type. The effective theory at the scale g 2 T is generated by a Boltzmann 
equation which includes the collision term for colour relaxation. 



1 Introduction 

At high temperature, the non-Abelian gauge theories describe weakly coupled 
plasmas, which, in a fust approximation, are very much alike the ordinary, elec- 
tromagnetic, plasmasEra. The plasma constituents, e.g., quarks and gluons for 
hot QCD, have typical momenta k ~ T, and take part in collective excitations 
which typically develop on a space-time scale A ~ (T is the tempera- 

ture, and g is the gauge coupling, assumed to be small.) Such excitations are 
similar to the familiar charge oscillations of the electromagnetic plasmasBLand 
can indeed be described by simple kinetic equations of the Vlasov type aid. 

But this simple analogy breaks down at the softer scale g 2 T. There, the 
non-Abelian plasmas enter a new non-perturbative regime where the coupling 
constant is small but the field strengths are large, so that perturbation theory 
breaks down because of large non-linear effects. Indeed, gluons obey Bosc- 
Einstein statistics, so the population of the soft (k <C T) gluon modes is 
strongly enhanced in thermal equilibrium: 

N a (k) ee -^—^ ~ - > 1, for k « T. (1) 

Thus, the long wavelength thermal fluctuations with A ~ l/g 2 T involve many 
quanta Nq ~ l/g 2 and behave in many respects as classical colour fields A£ 
with large amplitudes. We shall verify later that, typically, |A| = */ (A 2 ) ~ gT. 
This is a large fluctuation in the sense that the two terms in the soft covariant 
derivative D x = d x + igA are of the same order in g, d x ~ gA ~ g 2 T, so that 
the non-linear effects are indeed non-perturbative. In perturbation theory, such 



effects show up as infrared divergences in relation with the mutual interactions 
of the soft (k ~ g 2 T) magnetic gluons (cf. Sec. 2. a below). 

Moreover, infrared divergences are also associated with interactions among 
the hard (k ~ T) particles, as mediated by the exchange of soft magnetic 
gluons (or photons, in QED). This occurs, for instance, in the calculation of 
the quasiparticles damping rates in both QCD and QED (cf. Sec. 2.b). 

In order to deal with such problems, one has to go beyond ordinary pertur- 
bation theory. For an Abelian plasma, it is possible to eliminate the infrared 
problem of the damping rate by a specific resummation of the perturbation 
theory, based on the Bloch-Nordsieck (or eikonal) approximation u. For non- 
Abelian plasmas, the non-linear effects in the soft magnetic sector must be 
treated exactly, which requires numerical methods. An useful strategy in this 
sense — which is especially well suited for real-time calculations — is to first 
construct an effective theory for the soft fields, by integrating out the hard fields 
in perturbation theory, and then stud^ r J f bp 1 |effectiY/ : ' theory non-perturbatively, 
e.g., as a classical theory on a lattice aQi3'BE30'E3EjO. (The classical approx- 
imation will be discussed in Sec. 3 below.) This strategy can be seen as an 
extension to real time of the "dimensional reduction" generally performed in 
static calculations E£l. 

What I would like to show you here, is that kinetic theory is a powerful 
tool for constructing such effective theories. This has been first demonstrated 
for the collective dynamics at the scale gT, where we have shownta that simple, 
collisionless, kinetic equations resum an infinite number of one-loop diagrams 
with soft-external lines and hard loop momenta, the so-called "hard thermal 
loops" Ej'tZl. These equations will be reviewed below, in Sec. 4, together with 
the effective classical theory they generate 113, in Sec. 5. Then, in Sec. 6, I 
shall extend this approach to describe collective colour excitations at the softer 
scale g 2 T. Thisinvolves a Boltzmaim. equationEj-H which generates Bodeker's 
effective theory liil (see also Refs. E3o). 

2 Infrared problems in perturbation theory 

(a) Higher-order corrections to the free energy 

The infrared (IR) complications are most easily seen in the calculation of 
static quantities, like the pressure, as performed in the imaginary time for- 
malism^]. There, the energies of the gluon modes are purely imaginary and 

b I consider a purely Yang-Mills plasma; indeed, quarks are not important for the infrared 
physics to be discussed here. 



Figure 1: An n-loop ladder diagram contributing to the free energy in hot QCD. 



discrete, ko = iuj n = i2irnT, with integer n ( "Matsubara frequencies"), which 
provides a large "screening mass" \u> n \ ~ T for all the non-static (n ^ 0) 
modes. Thus, the IR problems can be associated only with the static (u> n = 0) 
Matsubara modes, and the most severe IR divergences are expected to come 
from diagrams involving the static modes alone. 

For instance, when computing higher order corrections to the free energy in 
hot QCD, one finds strong IR divergences from the ladder diagrams depicted in 
Fig. [l| In this diagram, all the propagators are static, and the loop integrations 
are three-dimensional. By power counting, this can be estimated ascl ~ 
<? 6 T 4 (g 2 T/fi) n (for n > 4 loops), where \i is an ad-hoc IR cutoff. 

For the electric gluons, an IR cutoff ~ gT is indeed generated dynamically, 
via Debye screening. In the magnetic sector, however, screening can only occur 
non-perturbatively, at the scale g 2 T. With /j, ~ g 2 T, all the ladder diagrams 
with four or more loops will contribute to the same order in g, namely, to order 
g 6 . Thus, perturbation theory breaks down, in the sense that we lose the usual 
connection between powers of the coupling constant and the number of loops. 

As mentioned in the Introduction, this breakdown is associated with long 
wavelength (A ~ l/g 2 T) thermal fluctuations with large amplitudes, A ~ gT. 
To see this, consider the free propagator of the magnetic gluon, in imaginary 
time: 

(A(r,x)A(0)) = T]T / " e-^^vj— J-. (2) 

„J Or k 2 + ui 

By letting r — > 0, x — > 0, and keeping only the contribution of the static modes 
(w n = 0) with soft momenta k ~ g 2 T, one obtains: 



(A 2 ) ^Tj^~ g *T 2 , 



(3) 



so that \A\ = (A 2 ) ~ gT, as anticipated. 

A different perspective on these IR problems follows by noticing that the 
diagram in Fig. [l] is actually a graph of three-dimensional QCD with coupling 
constant <?3 = gy/T. That is, the leading IR behaviour of hot QCD can be 
studied with the replacement (with A l a (x.) = T AtA\{t, x)) : 

Z 4 = J ZW£(r,x) rap j- jfV J d 3 x I F^F^ 

— Z 3 ^ | OAl(x) exp{-/3 1 d 3 xi^^j (4) 

which reduces the initial finite-temperature problem in four dimensions to an 
effective zero-temperature problem in three Euclidean dimensions. This is 
the crudest example of "dimensional reduction" , a strategy which consists 
in integrating out the non-static modes in perturbation theory in-order to 
obtain an effective three-dimensional theory for the static modes 113. This 
is convenient since three-dimensional lattice simulations are much easier to 
perform than in four dimensions. For instance, the magnetic mass, and the 
non-perturbative correction, of 0(g 6 ), to the pressure, have been computed in 
this wayE3. Moreover, the reduction to three dimensions has also pppiitted 
for some analytic non-perturbative studies of the magnetic screening til. 

(b) Quasiparticle damping rates 

Non-perturbative aspects, as related to the unscreened magnetic gluons, 
appear also in the computation of some real-time correlation functions. The 
most celebrated example in thai-sepse is the rate for anomalous baryon number 
violation at high temperaturel3olI2rlLij, or "hot sphaleron rate" , a quantity which 
has been extensively discussed at this conference (see, e.g., the contributions 
by P. Arnold, D. Bodeker, A. Krasnitz, G. Moore and L. Yaffe). 

Here, however, I would like to discuss an example which is conceptu- 
ally simpler (since it involves only the computation of a 2-point function), 
namely the calculation of the lifetime of the quasiparticlcs. A single particle 
excitation is created by adding a particle with momentum p (e.g., a trans- 
verse gluon) to the plasma initially in equilibrium. The added particle will 
then scatter off the other particles in the thermal bath (see Fig. |J), thus 
changing its momentum and colour. That is, the initial excitation will de- 
cay, and this can be measured from the corresponding retarded propagator 
Dn(t, p) = i6(i)([A(t, p), A(0, — p)]). A usual expectation is that D^(t, p) de- 
cays exponentially in time, \Dn(t, p)| 2 ~ e~ 2jt , which identifies the lifetime 
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Figure 2: Elastic scattering in the (resummed) Born approximation. The continuous lines 
refer to hard gluons, while the wavy line is the soft gluon exchanged in the collision. The 
blob indicates the resummation of the screening effects. 




of the single particle excitation as r = I/27. It turns ouLhpwever, that the 
calculation of 7 is afflicted with IR problems. SpecificallyE^'El : 

vMa,q)\ 2 + r^(go,g)l 2 |,(5) 

where *T>i(q) and *T> t (q) denote the propagators of the exchanged gluon in 
the electric and the magnetic channels, respectively. These are dressed so 
as to include the screening effects at the scale gT; indeed, with a bare gluon 
propagator, the integral in eq. (||) would be quadratically IR divergent, showing 
that the damping rate is dominated by soft momentum .transfers, q gT. The 
dressed propagators have the following IR behaviour Id-El (below, rtiu is the 
Debye mass, m 2 D = g 2 N c T 2 /3): 



*Vi(q - 0,«) * 1 a , *V t (q « q) ~ - . 1 2 , (6) 



which exhibits Debye screening in the electric sector and dynamical screen- 
ing in the magnetic sector. The latter is due to Landau damping El, i.e., the 
thermal absorbtion of the space-like gluons, which gives an imaginary part to 
the self-energy of the magnetic gluon: lm.Ht{qo -~{-Kq /Aq)m 2 D . This 

is proportional to the frequency qo since only the time-dependent (qo =/= 0) 
magnetic fields can transfer energy to the plasma constituents (as mechanical 
work), and thus get damped. 

Because of Debye screening, the electric contribution to the damping rate 
7/ is finite and of order g T 3 /m 2 D = 0(g 2 T). However, even after including 
the dynamical screening, the magnetic contribution 74 remains logarithmically 



divergent (below, /i is an ad-hoc IR cutoff) 



9 2 N C T [ m °dq g 2 N c T m D 

It - — : / — = — : In • (7) 

4?r q Ait fi 

The remaining divergence is associated to the static magnetic interactions, 
which are not screened at the scale gT. Assuming magnetic screening at the 
scale g 2 T, it follows that 7 ~ aN c T\n(l/g) (with a = g 2 /4n), where, however, 
the constant term under the logarithm cannot be determined (since sensitive to 
the non-perturbative screening mechanism). Thus, in QCD, the lifetime of the 
quasiparticles cannot be computed in perturbation theory beyond logarithmic 
accuracy. 

The same IR problem occurs also in QED, in the calculation of the lifetime 
of the charged particles. There, this is even, more intriguing, since thereJs no 
magnetic screening in the Abelian theories^. However, as shown in Ref. the 
Abelian problem can be solved by a further resummation of the perturbation 
theory, with the peculiar result that the electron propagator has an anomalous, 
non-exponential, decay law: Gn(t, p) ~ exp{— aTt ln(m_oi)} (see also Ref. £3). 



3 The classical approximation 

I come now to the main question to be addressed in this talk, which is, how to 
compute non-perturbative real-time correlations in hot gauge theories. Clearly, 
the standard lattice calculations, as formulated in imaginary time, are not ap- 
propriate for this problem. Fortunately, a fully quantum calculation is actually 
not needed: because of the Bose enhancement, the non-perturbative modes 
with k ~ g 2 T have large thermal occupation numbers (cf. eq. ([!])), which, 
by the correspondence principle, is a classical limit. For instance, the average 
energy per soft mode in thermal equilibrium, namely, 

elk) = ~ T for fc«T. (8) 

e p — 1 

is the same as expected from the classical equipartition theorem. Based on 
this observation, it has been suggested to compute the hot baryon numher 
violation rate through lattice simulations of a classical thermal field theoryErLj. 
The only question is, what is the correct classical theory ? 

It is well known that the classical approximation becomes meaningless 
at high momenta k <; T, where eq. (^) is not correct anylonger. At a first 
sight, one could expect this to be irrelevant for the problem at hand. Indeed, 
we are interested here in non-perturbative correlations, as determined by the 
dynamics at the scale g 2 T. Thus, one may expect such correlations not to be 



sensitive to the hard plasma modest! If this was the case, such quantities could 
be simply computed from the classical Yang-Mills theory at finite temperature, 
without worrying too much about its bad ultraviolet behaviour. 

But the previous examples in Sec. 2 show that this argument is too naive! 
: in the plasma, the soft and hard modes are coupled by the interactions, which 
results in screening effects which considerably modify the dynamics of the soft 
modes. For instance, eq. @ for *T> t (q) shows that, for large enough frequencies 
<7o, the soft magnetic fields are efficiently screened by Landau damping, and 
therefore decouple from the non-perturbative IR physics. For q ~ g 2 T, only 
the modes with very low frequencies qa £ q^ /m 2 D ~ g 4 T can take part in non- 
perturbative phenomena. According to Arnold, Son and YaffeB, this seta -the 
time scale for non-perturbative phenomena to be l/g 4 T (see also Refs. Erlldllj). 

The classical Yang-Mills theory does not describe correctly the screening 
effects due tft-ifte hard particles. For instance, it yields a linearly divergent 
Debye mass&Bll3 m 2 cl ~ g 2 TA, rather than the correct, quantum, result m 2 D ~ 
g 2 T 2 . A possible solution to this problem is to treat hard and soft modes on 
a different footing: first, the hard modes are integrated out in perturbation 
theory, which properly generates the screening corrections; then, the resulting 
effective theory for the soft modes is treated as a classical field theory, via non- 
perturbative methods (e.g., via lattice simulations). This strategy requires 
an unambiguous separation between hard and soft degrees of freedom, e.g., an 
intermediate cutoff fi, which moreover must be consistent with gauge symmetry 
and with the lattice implementation. If \i is chosen such as gT <&[ Ji,.<gLT, then 
one obtains the "hard thermal loop" (HTL) effective theory ^EME], to be 
presented in the next two sections. If, onJiiLrilocr hand, g 2 T <C [i -C gT, 
then one obtains Bodeker's effective theoryliJ'ESEJIla, to be discussed in Sec. 6 
below. 



4 Effective theories from kinetic equations 

To construct an effective theory for the soft modes, one needs to study the 
dynamics of the hard plasma constituents — here, transverse gluons with mo- 
menta k ~ T — in the presence of soft (q 5s gT) background fields A%. This is 
an off-equilibrium situation: the plasma is perturbed away from equilibrium by 
the background fields which induce long wavelength (A <; 1/gT) fluctuations 
in the colour density of the hard particles. Since ) » f (where f ~ 1/T is 
the mean interparticle distance), these are collective colour excitations. Since, 
furthermore, A S> At (where At = 1/fc ~ 1/T is the thermal wavelength of the 
rd particles), we expect such excitations to be described by kinetic theory 
And, indeed, the equations to be presented below □ can be viewed as a 



generalization of the Vlasov equation for ordinary plasmas 0. 

Specifically, the longwavelength colour excitations of the hard transverse 
gluons are described by a colour density matrix iV a j,(k, x) which, to the order 
of interest, can be written in the form: 

iV ab (k, x) = N (e k )S ab - gWabix, v) (dN Q /de k ), (9) 

where No(ek) = l/(e' 3£fc — 1) is the equilibrium distribution (with £/. = |k|), 
and the function W(x, v), which parametrizes the off-equilibrium deviation, is 
a colour matrix in the adjoint representation, W(x,~v) = W a (x,~v)T a , which 
depends upon the velocity v = k/e^ (a unit vector), but not upon the magni- 
tude k = |k| of the momentum. It satisfies the following simple equation: 

(v-D x ) ab W b (x,v) = v^x), (10) 

where = (1, v), = <9 M + igA%T a is the covariant derivative defined by 
the background field, and E z a is the chromoelectric mean field. The system is 
closed by the Yang- Mills equations for the soft fields A%, namely: 

(D v F^) a (x) = j£(x), (11) 

with the induced current: 

j^x)=m 2 D J^v^W a (x,v), (12) 

where the angular integral J dSl runs over the orientations of v. 

Eq. ([To]) is a collisionless kinetic equation. It has been obtaineda from the 
general Dyson-Schwinger equations for the off-equilibrium plasma, by neglect- 
ing the collisions among the hard particles and by performing a gauge-covariant 
gradient expansion which takes profit of the assumed separation of scales. Note 
that all these approximations are controlled by the same small parameter, the 
coupling strength g, so that eq. ( |l0| ) is actually correct to leading order in g. 

By formally solving eq. (|To|), we can express the current in terms of the 
gauge fields and thus obtain an effective Yang-Mills equation which in- 
volves the soft fields alone: 

Eq. (|l^) describes the propagation of soft colour fields in the high-T plasma. 
The hard particles are not explicit anymore, since they have been integrated to 
yield the induced current in the r.h.s. By expanding this current in powers of 



the gauge fields one generates^ all the HTL's of Braaten and PisarskiB, which 
encompass the screening phenomena at the scale gT (cf. Sec. 2). However, 
because of the non-local structure of the current (note the covariant derivative 
in the denominator), eq. (13) is not very convenient for the construction of the 
classical thermal theory, to which I turn now. 

5 The classical effective theory 



I shalL-uow use eqs. (|HJ)-(|12|) to define a classical field theory at finite temper- 
atureEU. As usual with gauge theories, this is most easily done in the temporal 
gauge Aq — 0, where the equations read: 

Ef = -d Q A?, 

-d E? + e ijk (D j B k ) a =m 2 D J^ Vl W a (x, v), 

(8q + v • D) ofc Wf, — v • E a , (14) 

together with the constraint expressing Gauss' law (i.e., the /i = component 
ofeq. dip) : 

G a (x) = (D • E) a + m 2 D [ ^ W a (x, v) = 0. (15) 

J 4tt 



Note that eqs. (|14|) are not in canonical form: this is already obvious from the 
fact that we have an odd number of equations. Still, it can be verified that these 
equations are conservative; thejeorresponding, conserved energy functional has 
the gauge-invariant expression au : 

H= \ [ d 3 4E a -E Q + B a -B a + m 2 D [ ^ W a (x, v) W a {x, v)l. (16) 



2 



Eqs. (|l4T)-([l6l) define an effective theory for the soft degrees of freedom. Besides 
the soft colour fields A z a (x) and E l a (x), these equations also involve the auxiliary 
fields W a (x,v) which simulate the hard thermal gluons (or, more precisely, 
their long wavelength colour fluctuations). We are interested in computing the 
(non-perturbative) real-time correlations of the fields A l a . To this aim, we need 
to construct the thermal partition function for this classical field theory. 

Recall first how this is done in some generic theory: the thermal expecta- 
tion values are obtained by first solving the classical equations of motion for 
given initial conditions, and then averaging over the initial conditions (i.e., over 
the classical "phase-space") with the Boltzmann weight exp(—(3H). Since the 
initial conditions (say <^(x) and </>(x) for a scalar theory) depend only on the 



spatial coordinate x, the phase space integration is actually a three-dimensional 
functional integral, which can be implemented on a lattice in the standard way. 

For the problem at hand, the classical phase-space is determined by the 
initial conditions to eqs. (|l4|), that is, 

Al(0,x)=^(x), ££(0,x)=£*(x), W.(0 > x > v)=W o (x,v), (17) 

and the thermal weight is provided by the Hamiltonian in eq. (|l6|). Then, the 
(real-time) thermal correlation functions of the fields A l a can be obtained from 
the following generating functional: 

Z cl [J?] = J VE?VA a l VW a 5{g a ) exp j-/3H + j d 4 zJ? [x)A a l {x]^ , (18) 

where A* (t,x) is the solution to eqs. ( O ) with the initial conditions (p"7|), and 
Q a and 7i are expressed in terms of the initial fields, cf. eqs. ( |l5| ) and (fl6|). 
Physically, the fluctuations in the initial conditions W(x, v) for the auxiliary 
fields can be interpreted as a thermal noise due to the hard particles and which 
drives the soft fields toward thermal equilbrium. It is this noise which generates 
the Landau damping of the soft correlation functions E3. 

In particular, for J = 0, eq. ( |T8] ) yields the following expression for the 
free energy of the classical thermal radiation: 

Z cl = J VA a VA«exp {-/? J d 3 x Q(^-) 2 + \{ViA»f + ^(Ag) 2 ) } (19) 



which is also the result expected from dimensional reduction (see eq. 

To complete the construction of the effective theory, eq. (|T|) must be 
supplemented with an UV cutoff fj,, which is the scale separating hard from 
soft degrees of freedom: gT <C <C T. Correspondingly, the Hamiltonian 
must be extended to include /i-dependent counterterms, chosen so as to cancel 
the cutoff dependence of the classical theory in any complete calculation. In 
practice, however, this "matching" turns out to be difficult to achieve E3u, 
mainly because of the constraints of the lattice implementation. It is therefore 
rewarding that the results of the first lattice simulations lJ of eqs. (|l^-([l8"|) 
appear to be quite robust and insensitive to lattice artifacts. In particular, 
the hot sphalcron rate obtained in this way is consistent with the previous 
calculations in Ref. 



6 A Boltzmann equation for colour 

Since the "semi-hard" modes with q ~ gT are also perturbative, it is possible 
to integrate them out as well, and thus get an effective theory involving only 



the "ultrasoft" modes with q ~ g 2 T0. Then, the collisions among the hard 
particles, as mediated by the semi-hard fields, must be included explicitely in 
the kinetic equation. That is, the previous equation (|l^) must be generalized 
so as to include the collisions terms. This is also necessary for consistency: 
for colour fluctuations at the scale g 2 T, the effects of the collisions_among the 
plasma particles become as important as those of the mean fields E_3. 

The resulting kinetic equation is a Boltzmann equation describing the 
propagation and relaxation of longwavelength (A ~ 1 / g 2 T]_rplQm- excitations. 
It allows, in particular, to compute thexplour conductivitylI3E3'E3, thus clarify- 
ing some previous work on this subjectrJ. To leading logarithmic accuracy (see 
below), this equation has been first obtained by BodekerO. Then, Arnold, Son 
and Yaffe have shown 113 that Bodeker's theory can be generated by a rather 
simple Boltzmann equation—which has been physically motivated, but not 
rigorously-proven, in Refs. 113. Recently, we have given a derivation of this 
equationl!3, starting from the quantum field equations. This has clarified the 
nature of the approximations involved, thus fixing its range of applicability. 

Remarkably, even after the inclusion of the collision term, the density ma- 
trix 7V a b(k, x) preserves the same structure as in the mean field approximation 
(cf. eq. (Jsj) ) , but the functions W (x,v) satisfy a more complicated equation 
(compare to eq. ([lO])) : 

f ($Cv ■ v')W a (x v')) 1 

(v ■ D x y b W b (x, v) = v • E a (x) 7 [W*(x, v) - [ 1 ^ ■ ^ ,} j(.20) 

The new feature here is the collision term in the r.h.s. This is proportional 
to the quasiparticle damping rate 7 (cf. Sec. 2.b), which thus appears to set 
the scale for the colour relaxation time: r co ; ~ 1/7 ~ l/(g 2 Tln(l/g)). The 
other notations above are as follows: the angular brackets in the collision term 
denote angular average over the directions of the unit vector v' (as in eq. (22) 
below), and the quantity $(v ■ v') is given by: 



*( v ' O = / T^f %0 - q • v)S(q - q ■ v') *©,(«) + (v t ■ v' t ) *V t (q) .(21) 



2 

- n • vML - n-v'l *T),(n\ 4- (v. ■ v'A *T)An\ 

with *T>i(q) and *T> t (q) defined after eq. (||). Up to a normalization, $(v • v') 
is the total interaction rate for two hard particles with momenta k and p 
(and velocities v = k and v' = p) in the (resummed) Born approximation, as 
illustrated in Fig. ^ (v t and are the transverse projections of the velocities 
with respect to the momentum q of the exchanged gluon). The damping rate 
7 is obtained from <i>(v • v') as follows (cf. eq. (|t])) : 

t = —t~ J $(v ' v " H1/9) - (22) 



There are two important remarks about the previous equations: 

First, whereas most transport phenomena are dominated by large m£t: 
mentum transfers gT ^ q ^ T, so that the typical relaxation times are E3 
r tr ~ 1 / (g A T / g)) , the relaxation of colour excitations turns out to be 
dominated by soft gluon exchanges, g 2 T ^ q ^ gT,^s_tha quasiparticle damp- 
ing rate. There is a simple physical reason for thatEjOO : unlike momentum 
fluctuations, which require a large angle scattering to relax, colour can be ef- 
ficiently exchanged .in. airy: scattering, even a small angle one. This yields a 



colour conductivity tallinla er c ~ T/ln(l/g), to be contrasted with the usual, 
electric, conductivity E3 : a e i ~ T/(e 2 ln(l/e)). Note also that it is the same 
physical process — namely, the scattering via one-gluon exchange in Fig. — 
which provides relaxation for both single-particle and collective (momentum 
or colour) excitations. If, nevertheless, the relevant time scales turn out not 
to be the same (namely r ~ t~j 3> T tr ), it is because of specific cancellations 
among various collision terms c3, which occur in the calculation of most trans- 
port coefficients, but not in that of the quasiparticle lifetime r, or in that of 
the relaxation time of colour excitations r co /. 

Second, strictly speaking, the functional form of the collision term in 
eq. ( |2(i| ) is only valid to leading logarithmic accuracy, because of the approx- 
imations performed in its derivation. This limitation comes from the poor 
convergence of the gradient expansion when the range of the interactions be- 
comes comparable to the scale of the systemJxihomogeneitiesEa. Within this 
logarithmic accuracy, eq. ( |20| ) can be shownE3t3 to reduce to Bodeker's equa- 
tionta. 

Note finally that, even though conceptually interesting, the leading loga- 
rithmic approximation appears to be of little use for the calculation jOf-the hot 
baryon number violation rate. Indeed, the numerical calculations L2t2c2 show 
that, for realistic values of g, the constant term under the logarithm is sensibly 
larger than \n(l/g). 
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